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ABSTRACT 

We present analyses of a 50 ks observation of the supergiant X-ray binary system CygnusX-l/HDE 226868 
taken with the Chandra High Energy Transmission Grating Spectrometer (HETGS). CygX-1 was in its spec- 
trally hard state and the observation was performed during superior conjunction of the black hole, allowing for 
the spectroscopic analysis of the accreted stellar wind along the line of sight. A significant part of the observation 
covers X-ray dips as commonly observed for CygX-1 at this orbital phase, however, here we analyze only the 
high count rate nondip spectrum. The full 0.5-10keV continuum can be described by a single model consisting 
of a disk, a narrow and a relativistically broadened Fe Ka line, and a power-law component, which is consistent 
with simultaneous RXTE broad band data. We detect absorption edges from overabundant neutral O, Ne, and Fe, 
and absorption line series from highly ionized ions and infer column densities and Doppler shifts. With emission 
lines of He-like Mg XI, we detect two plasma components with velocities and densities consistent with the base 
of the spherical wind and a focused wind. A simple simulation of the photoionization zone suggests that large 
parts of the spherical wind outside of the focused stream are completely ionized, which is consistent with the low 
velocities (<200kms -1 ) observed in the absorption lines, as the position of absorbers in a spherical wind at low 
projected velocity is well constrained. Our observations provide input for models that couple the wind activity of 
HDE 226868 to the properties of the accretion flow onto the black hole. 

Key words: accretion, accretion disks - stars: individual (HDE226868, CygX-1) - stars: winds, outflows - 
techniques: spectroscopic - X-rays: binaries 



1. INTRODUCTION 

Cygnus X-l was discovered in 1964 (Bowyer et al. 1965) and 
soon identified as a high-mass X-ray binary system (HMXB) 
with an orbital period of 5.6 d (Murdin & Webster 1971; 
Webster & Murdin 1972; Bolton 1972). It consists of the su- 
pergiant 09.7 star HDE 226868 (Walborn 1973; Humphreys 
1978) and a compact object, which is dynamically constrained 
to be a black hole (Gies & Bolton 1982). The detailed spectro- 
scopic analysis of HDE 226868 by Herreroetal. (1995) gives 
a stellar mass M* w 18M Q , leading to a mass of M BH ^1OM 
for the black hole, if an inclination i sa 35° is assumed. Note 
that Ziolkowski (2005) derives a mass of M* = (40 ± 5)M Q 
from the evolutionary state of HDE 226868, corresponding to 
M BH = (2O±5)M , while Shaposhnikov & Titarchuk (2007) 
claim Mbh = (8.7±0.8)M Q from X-ray spectral-timing rela- 
tions. 

Cyg X-l is usually found in one of the two states that are dis- 
tinguished by the soft X-ray luminosity and spectral shape, the 
timing properties, and the radio flux (see, e.g., Pottschmidt et al. 
2003; Gleissner et al. 2004a,b; Wilms et al. 2006): the low/hard 
state is characterized by a lower luminosity below lOkeV, a 
hard Comptonization power-law spectrum (photon index T ~ 
1.7) with a cutoff at high energies (folding energy Zsfdd ~ 
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150keV) and strong variability of ^30% root mean square 
(rms). Radio emission is detected at the ~15mJy level. In 
the high/soft state, the soft X-ray spectrum is dominated by a 
bright and much less variable (only few % rms) thermal disk 
component, and the source is invisible in the radio. Within the 
classification of Remillard & McClintock (2006), the high/soft 
state of Cyg X-l corresponds to the steep power-law state rather 
than to the thermal state, as a power-law spectrum with pho- 
ton index T ~ 2.5 may extend up to ^-TOMeV (Zhang et al. 
1997; McConnell et al. 2002; Cadolle Bel et al. 2006). Most 
of the time, Cyg X-l is found in the hard state, but transitions 
to the soft state and back after a few weeks or months are 
common every few years. Transitional or intermediate states 
(Belloni et al. 1996) are often accompanied by radio and/or X- 
ray flares. Similar to a transition to the soft state, the spec- 
trum softens during these flares and the variability is reduced. 
This behavior is called a "failed state transition" if the true soft 
state is not reached (Pottschmidt et al. 2000, 2003). Transi- 
tional states have occurred more frequently since mid- 1999 than 
before (Wilms et al. 2006), which might indicate changes in the 
mass-accretion rate due to a slight expansion of HDE 226868 
(Karitskaya et al. 2006). 

HMXBs are believed to be powered by accretion from the 
stellar wind. The accretion rate and therefore X-ray luminos- 
ity and spectral state are thus very sensitive to the wind's de- 
tailed properties such as velocity, density, and ionization. For 
HDE 226868, Gies et al. (2003) found an anticorrelation be- 
tween the Ha equivalent width (an indicator for the wind mass 
loss rate M*) and the X-ray flux. Considering the photoioniza- 
tion of the wind would allow for a self-consistent explanation 
(see, e.g., Blondin 1994): a lower mass loss gives a lower wind 
density and therefore higher degree of ionization due to the 
irradiation of hard X-rays, which reduces the driving force of 
HDE 226868 's UV photons on the wind and results in a lower 
wind velocity v, leading finally to a higher accretion rate (oc 
M*/v 4 , Bondi & Hoyle 1944). However, Gies et al. (2008) find 
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suggestions that the photoionization and velocity of the wind 
might be similar during both hard and soft states. UV obser- 
vations allow the photoionization in the HDE 226868 /Cyg X-l 
system to be probed: Vrtilek et al. (2008) reported P Cygni pro- 
files of N V, C IV, and Si IV with weaker absorption compo- 
nents at orbital phase <fi OI b w 0.5, i.e., when the black hole 
is in the foreground of the supergiant. This reduced absorp- 
tion, which was already found by Treves et al. (1980), is due to 
the Hatchett & McCray (1977) effect, showing that those ions 
become superionized by the X-ray source. Gies et al. (2008) 
model the orbital variations of the UV lines assuming that the 
wind of HDE 226868 is restricted to the shadow wind from 
the shielded side of the stellar surface (Blondin 1994), i.e., the 
Stromgren (1939) zone of Cyg X-l extends to the donor star. 
However, this assumption applies only to the spherical part of 
the wind, which might therefore hardly contribute to the mass 
accretion of Cyg X- 1 . 

As HDE 226868 is close to filling its Roche lobe (Conti 1978; 
Gies & Bolton 1986a,b), the wind is not spherically symmetri- 
cal as for isolated stars, but strongly enhanced toward the black 
hole ("focused wind"; Friend & Castor 1982). The strongest 
wind absorption lines in the optical are therefore observed at the 
conjunction phases (Gies et al. 2003). Similarly, X-ray absorp- 
tion dips occur preferentially around <f> OI b = 0, i.e., during su- 
perior conjunction of the black hole (Balucihska-Church et al. 
2000). These dips are probably caused by dense, neutral 
clumps, formed in the focused wind where the photoioniza- 
tion is reduced, although recent analyses have also suggested 
that part of the dipping activity may result from the interac- 
tion of the focused wind with the edge of the accretion disk 
(Poutanen et al. 2008). 

The photoionization and dynamics of both the spherical and 
focused winds can also be investigated with the high-resolution 
grating spectrometers of the modern X-ray observatories Chan- 
dra or XMM-Newton. As none of the previously reported ob- 
servations of Cyg X-l was performed at orbital phase <p 0I b = 
and in the hard state, when the wind is probably denser and less 
ionized than in the soft state, the Chandra observation presented 
here allows for the most detailed investigation of the focused 
wind to date. 

The remainder of this paper is organized as follows: in Sec- 
tion 2, we describe our observations of Cyg X-l with Chandra 
and the Rossi X-Ray Timing Explorer (RXTE), and how we 
model CCD pile-up for the Oiandra-HETGS data. We present 
our investigations in Section 3: after investigating the light 
curves, we model the nondip continuum and analyze neutral 
absorption edges and absorption lines from the highly ionized 
stellar wind - and the few emission lines from He-like ions, 
which indicate two plasma components. In Section 4, we dis- 
cuss models for the stellar wind and the photoionization zone. 
We summarize our results after comparing them with those of 
the previous Chandra observations of Cyg X-l. 

2. OBSERVATION AND DATA REDUCTION 

2.1. Chandra ACIS-S/HETGS Observation 

Cyg X-l was observed on 2003 April 19 and 20 by the Chan- 
dra X-Ray Observatory, see Table 1 . An overview on all its in- 
struments is given by the Proposers' Observatory Guide (CXC 
2006). In the first four months of 2003, the RXTE All-Sky 
Monitor (ASM; see Doty 1994; Levine et al. 1996) showed the 
source's 1.5-12keV count rate to be generally below 50 ASM- 
cps (Figure 1). At the time of our Chandra observation, it was 
less than 25cps (0.33 Crab), typically indicative of the source 
being in its low/hard state (Wilms et al. 2006). 



Table 1 

Observations of Cyg X-l 



Satellite / 


Start 


Stop 


Exposure" 


Count Rate" 


Instrument 


(MJD) 


(MJD) 


(ks) 


(cps) 


Chandra 


52748.70 


52749.28 


(47.2) 


( ■■■ ) 


MEGil* 


(52748.70) 


(52749.28) 


16.1 


2x27 


HEG±1 C 


(52748.70) 


(52749.28) 


16.1 


2 X 17 


RXTE 


52748.08 


52749.18 






PCA rf 


52748.74 


52748.78 


3.0 


1456 


HEXTE a+b c 


52748.74 


52748.78 


1.1 


2x186 



Notes. 

" For the Chandra data, the nondip GTIs (see Figure 4) have been used. 

For RXTE, the 1 1th orbit was considered. 
* The Oiandra-MEG spectra cover pa 0.8-6 keV (1-99 % quantiles). 
c The C/iandra-HEG spectra cover fa 1-7.5 keV (1-99 % quantiles). 
d The RXTE-PCA data from 4 to 20 keV has been used. 
e The RXTE -HEXTE data from 20 to 250 keV has been used. 
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Figure 1. Brightness of Cyg X-l as seen by the ASM on board RXTE. During 
the Chandra observation (marked by a line), the source was still in its low/hard 
state: the 1.5-12keV count rate did not exceed 25 cps. In spite of the high 
intrinsic variability, the high/soft state during June, July, and August can clearly 
be distinguished, a result which is also found by spectral analysis (Wilms et al. 
2006). 

The High Energy Transmission Grating Spectrometer 
(HETGS), containing high and medium energy gratings 
(HEG/MEG; see Canizares et al. 2005) was used to disperse X- 
ray spectra with the highest resolution (CXC 2006, Table 8.1): 

AAheg = 5.5mA and AA M eg = 1 1 mA (1) 

As only half of the spectroscopy array of the Advanced 
CCD Imaging Spectrometer (ACIS-S; see Garmire et al. 2003), 
namely a 512 pixel broad subarray, was operated in the timed 
event (TE) mode, the six CCDs could be read out after expo- 
sure times of ?f rame = 1 .7 s and the position of each event is well 
determined. Photons from the different gratings can thus easily 
be distinguished due to the different dispersion directions of the 
HEG and the MEG. Even in its low state, however, Cyg X-l is 
so bright that several photons may pile up in a CCD pixel during 
one readout frame. Both events cannot be discriminated and are 
interpreted as a single photon with larger energy. As the undis- 
persed image would have been completely piled up, only 10% 
of those events in a 40 x 38 pixel window have been transmit- 
ted in order to save telemetry capacity. The first-order spectra 
are, however, only moderately affected, which can be modeled 
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very well (Section 2.2). The alternative to this approach would 
have been to use the continuous clocking (CC) mode, where 
the ACIS chips are read out continuously in 2.85 ms, but only 
the position perpendicular to the readout direction can be de- 
termined for every photon event. The CC mode was, however, 
avoided due to difficulties in the reconstruction of HEG and 
MEG spectra and other calibration issues. 

The undispersed position of the source is required for the 
wavelength calibration of the spectra. We redetermined it to 
R.A. = 19 h 58 m 21. s 67, 8 = +35° 12'5"83 from the intersection 
of the HEG and MEG arm and the readout streak (Ishibashi 
2006). Afterward, the event lists were reduced using the stan- 
dard software from the Chandra X-ray Center (CXC), CIAO 3.3 
with CALDB3.2.3. 8 Exceptionally narrow extraction regions 
had to be chosen as the background spectrum would other- 
wise have been dominated by the dispersed extended X-ray- 
scattering halo around the source (Xiang et al. 2005). The fur- 
ther analysis was performed with the Interactive Spectral Inter- 
pretation System (ISIS) 1.4.9 (Houck 2002). 9 

We use the four first-order MEG and HEG spectra (with two 
dispersion directions each, called +1 and -1 in the following) 
which provide the best signal-to-noise ratio (S/N). The "second- 
order spectra" are dominantly formed by piled first-order events 
which reach the other order sorting window of data extraction 
(defined in dispersion-energy space) when the energy of two 
first-order photons accumulates. This effect is most evident for 
the MEG, whose even dispersion orders are suppressed by con- 
struction of the grating bars (CXC 2006). 

2.2. Model for Pile-Up in Grating Observations 

For the first-order spectra, pile-up causes a pure reduction of 
count rate: a multiple event, i.e., the detection of more than 
one photon in a CCD pixel during one readout time, which 
cannot be separated, is either rejected by grade selection dur- 
ing the data processing or migrates to a higher-order spectrum. 
The Poisson probability for single events in a 3 x 3 pixel event- 
detection cell i (see Davis 2002, 2003; CXC 2005) is 

Pdi) = A(i)-exp(-A(i)) , (2) 

where the expected number of events, A(z) = 70 • C tot (i), is given 
by the total spectral count rate, C to t0X at this position (in units 
of counts per A and s), and where the constant 70 is 

70 = 3 AA • ff rame , (3) 

where A A is the resolution of the spectrometer of Equation (1), 
and ff rame is the frame time. We therefore describe the pile-up 
in the first-order spectra with the nonlinear convolution model 
simple_gpile2 in ISIS, 10 which exponentially reduces the 
predicted count rate C(A) according to Equation (2): 

C'(A) = C(A)-exp(- 7 -C tot (A)) (4) 

Here, the scale 7 w 70 is left as fit parameter, and C to t(A) also 
takes the photons into account which are dispersed in a higher 
order m < 3. The count rates are estimated from the correspond- 
ing effective areas A,„ and the assumed photon flux S: 

C,ot(A) = V 3 A m (X/m) ■ S(X/m) (5) 

simple_gpile2 is based on the simple_gpile model 
(CXC 2005; Nowak et al. 2008), which parameterizes the 

8 See http://cxc.harvard.edu/ciao3. 3/. 

9 See http://space.mit.edu/cxc/isis/. 

10 The ISIS/S-LANG code for simple_gpile2 is available online at 
http://pulsar.sternwarte.uni-erlangen.de/hanke/X-ray/code/simple_gpile2.sl. 
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Figure 2. Pile-up in the HEG (black curves) and MEG (gray curves) data: 
the lower spectra, which are shifted by a factor of 0.4, show the uncorrected 
data. The dashed line shows the model (free of pile-up). The MEG spec- 
trum suffers from significant pile-up losses around 2keV, where the highest 
count rate is obtained. The upper spectra show the pile-up corrected spectra. 
Note that we show ISIS ' (model independent) flux spectra only for illustration; 
simple_gpile(2) operates on the count rates predicted by a model. 

strength of pile-up by the (maximum) pile-up fraction 
p = 1 -exp (-7 ■ max{C t ot}), while using the parameter 7 of 
simple_gpile2 avoids to have a nonlocal model which de- 
pends on the flux at the position of the highest pile-up. 

The effect of pile-up is stronger in the MEG spectra than in 
the HEG spectra due to the lower dispersion and higher effective 
area of the MEG. The apparent flux reduction is most significant 
near 2 keV where the spectrometer has the largest efficiency and 
the highest count rates are obtained (see Figure 2). It can, how- 
ever, excellently be modeled with simple_gpile2. When 
fitting a spectrum, there is always a strong correlation between 
the pile-up scale 7 and the corresponding flux normalization 
factor, e.g., the relative cross-calibration factor c introduced in 
Section 3.2. The best-fit values for the pile-up scales 7 found 
in our data analysis (see Table 2) are only slightly larger than 
expected from Equation (3) (namely 7o.meg±i = 5.6 x 10~ 2 s A 
and 7o.heg±i = 2.8 x 10~ 2 s A), and the calibration factors c are 
consistent with 1, except for the HEG-1 spectrum, for which 
both the largest 7 - 70 and c were found. Given the presence 
of the j—c correlation, we consider the latter to be a numerical 
artifact. 

According to the simple_gpile2 model, the MEG+1 
spectra suffer from >30 % pile-up for 6 A < A < 8.3 A, peaking 
at Pmeg+i = 45 % in the Si XIII f emission line at 6.74 A. Except 
for some emission lines, among them the Fe Kq line, the con- 
tinuum pile-up fraction of the HEG spectra is below 17 %. For 
the HEG+1 spectrum, the reduction is less than 10% outside 
the ranges 2.09 A < A < 4.08 A and 6.05 A < A < 6.93 A. 

2.3. RXTE Observation 

While Chandra's HETGS provides a high spectral resolution, 
the energy range covered is rather limited. Within the frame- 
work of our RXTE monitoring campaign, the broadband spec- 
trum of Cyg X-l was measured regularly, i.e., at least biweekly, 
since 1998 (see Wilms et al. 2006, and references therein). The 
observation on 2004 April 19 was extended to provide hard X- 
ray data simultaneously with the Chandra observation. More 
than one day was covered by 17 RXTE orbits of ~47 min good 
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Figure 3. Coverage of the simultaneous RXTE and Chandra observations. The plot shows the background-subtracted light curves with a time resolution of 64 s on a 
separate y axis each. Top: RXTE -HEXTE (20-250 keV). The count rates have been corrected for the detector dead time. Center: RXTE-PCA (4-20 keV), normalized 
by the number of active PCUs. Bottom: C/iandra-HETGS (0.5-12 keV, only first-order events; see Figure 4 for more details). The numbers label the RXTE orbits, 
vertical lines mark the first part of the nondip spectrum (see Section 3 . 1 and Figure 4), which completely covers RXTE orbit 1 1 . 

orbital phase of the binary system 
0.94 0.96 0.98 0.00 0.02 




20 30 
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Figure 4. Top: full 0.5-12 keV band Chandra light curve. Bottom: ratio of 
0.7-1 keV band and 2.1-7.2keV band count rates. Absorption dips - at first 
compact, then with complex substructure - show up with a reduced flux and 
spectral hardening. Count rates >82.7 cps define the nondip data (dark). 

time, interrupted by ~-49min intervals when CygX-1 was not 
observable due to Earth occultations or passages through the 
South Atlantic Anomaly (SAA), see Figure 3. 

Data in the 4-20 keV range from the Proportional Counter 
Array (PCA; Jahoda et al. 1996) and in the 20-250 keV range 
from the High Energy X-Ray Timing Experiment (HEXTE; 
Gruberet al. 1996) were used. The data were extracted using 
HEASOFT 6.3.1, 11 following standard data screening proce- 
dures as recommended by the RXTE Guest Observer Facility. 
Data were only used if taken more than 30 minutes away from 
the SAA. For the PCA, only data taken in the top layer of the 
proportional counter were included in the final spectrum, and 
no additional systematic error was added to the spectrum. Dur- 
ing the observation, different sets of Proportional Counter Units 
(PCUs) were operative. During the 1 1th orbit, extensively used 
in this work, PCUs 1 and 4 were off. 

3. ANALYSIS 

3.1. Light Curve 

The Chandra observation covers a phase range between 
~0.93 and ~0.03 in the 5.599829d binary orbit (Giesetal. 
2003, whose epoch is HJD 2451730.449 ± 0.008). The top 
panel of Figure 4 shows the light curve of first-order events 
(MEG±1, HEG±1) in the full band accessible with Chandra- 
HETGS. During several dip events, the flux is considerably 

11 See http://heasarc.gsfc.nasa.gov/lheasoft/. 
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Figure 5. RXTE 4-250 keV broadband continuum spectrum, as measured with 
PCA below 20 keV and above 20 keV with HEXTE, which can be described by 
a broken power-law with high-energy cutoff and a weak iron line, see Table 2. 
Due to the joint fit with Chandra, the Fe line consists of narrow and a broad com- 
ponent, see Figure 6. The HEXTE spectra shown are renormalized to match the 
PCA flux, as the absolute calibrations of these two instruments differ by ~ 15%. 

reduced. The absorption dips distinguish themselves also by 
spectral hardening (the bottom panel of Figure 4). We extract 
a 16. Iks nondip spectrum, which is the subject of this paper, 
from all times when the total count rate exceeds 82.7 cps; these 
are indicated by dark points in Figure 4. The analysis of dip 
spectra will be described in a subsequent paper. 

The RXTE light curve shows considerable variability as well. 
Although the dips are more obviously detected with Chandra 
in the soft X-ray band, similar structures are also seen with 
RXTE-PCA or even -HEXTE, especially in the last RXTE or- 
bits of these observations (Figure 3). We chose to infer the 
nondip broadband spectrum from the RXTE data taken during 
the 11th orbit, which was performed entirely during the first 
part of the nondip phase. Other parts are interrupted by dips, 
occultations, or have nonuniform PCU configuration. 

3.2. Continuum Spectrum 

The broadband spectrum of Cyg X-l in the hard state can be 
described by a broken power-law with exponential cutoff (Fig- 
ure 5). Since the parameters of this phenomenological model 
are correlated with those of physical Comptonization models 
(see, e.g., Wilms et al. 2006, Fig. 11), we are justified in using 
the aforementioned simple continuum model for this paper fo- 
cusing on the spectroscopy of the wind. We describe the whole 
0.5-250 keV spectrum consistently with one broken power-law 
model, i.e., there is no need to fit the continuum locally. 



CHANDRA X-RAY SPECTROSCOPY OF THE FOCUSED WIND IN THE CYG X-l SYSTEM. I. 



Table 2 

Fit Parameters of the Continuum Nondip Hard State Spectrum of Cyg X- 1 







Fit to the 


Joint Fit to Both the 


Parameter 


Unit 


Chandra 
Spectra Only 


Chandra and RXTE 
Spectra 


Photoabsorption 








iV H 


10 21 cm" 2 


3.52 ±0.04 


5.4 ±0.4° 


(Broken) power-law 








Fi (HETGS) 




1.51 ±0.01 


1.60±0.01 


r, (pc A) 






1.73±0.01 


^break 


keV 




9.0+0; 3 


r 2 






1.50±0.01 


norm 


s" 1 cm" 2 keV" 


1.15±0.01 


1.33 ±0.03 


High-energy cutoff 








Eeui 


keV 




24+| 


Efold 


keV 




204 ±9 


Disk black body 








/tdisk (Equation 6) 


10 3 




23+17 


kT co \ 


keV 




25+ 03 

u -0.02 


Narrow iron Ka line 








^O.narrow 


keV 




6 4001+ 00004 


& narrow 


eV(!) 




<0.1 


^narrow 


10" 3 s" 1 cm" 2 




1.0±0.3 


Broad iron Ka line 








^O.broad 


keV 




6.4 (fixed) 


""broad 


keV 




0.6±0.1 


Abroad 


10" 3 s" 1 cm" 2 




4 3+ 2 - 7 

+- J -0.5 


Relative flux calibration (constant factor) 




CMEG-l 




1 (fixed) 


1 (fixed) 


C MEG+1 




1.00 ±0.01 


0.99 ±0.01 


CHEG-l 




1.04 ±0.01 


1.04±0.01 


C HEG+1 




1.01 ±0.01 


1 AT I fl Al 

l .02 ± 0.0 1 


CPCA 






i i g -i- n m 

1 . 1 o zt \>A>3 


CHEXTE 






1 rw+0.02 


Pile-up scales 








7MEG-1 


10" 2 s A 


5.6±0.1 


5 6 ± 1 


TMEG+i 


10" 2 s A 


6.0±0.1 


5 9 ± 1 


7HEG-1 


10" 2 s A 


4.5 ±0.3 


4.6 ±0.3 


7HEG+1 


10" 2 s A 


3.5 ±0.3 


3.8 ± 0.3 


Fit-statistics 








x 2 




11745 


12180 


dof 




11274*-10 c 


11274 6 + 293' J -25 C 


-^red 




1.04 


1.06 


Absorbed flux (pile-up corrected) 






So.5-10keV 


photons s 1 cm" 


2 


1.4 


fb.5-10kcV 


10"' erg s" 1 cm" 


-2 


7.4 


So.5-250keV 


photons s" 1 cm" 


2 


1.8 


^0.5-250keV 


10" 9 erg s" 1 cm 


-2 


24 


Unabsorbed luminosity, assuming d 


= 2.5 kpc (Ninkov et al. 1987a) 


^0.5-lOkeV 


10 37 ergs" 1 


0.67 


0.78 


^0.5-250 kcV 


10 37 ergs" 1 




3.1 



Notes. Error bars indicate 90% confidence intervals for one interesting 
parameter. 

" In the joint fit, photoabsorption was described with tbnew model 
(Wilms et al. 2000; Juett et al. 2006b) and the best-fitting abundances (see Section 

* All data have been rebinned to contain > 50 counts bin" 1 . 

c The model contains actually more parameters, as the absorption lines have 
already been included (see Tables 4-6). 




Energy [keV] 

Figure 6. Oiandra-HEG spectrum in the Fe line region. The model includes 
both the narrow and the broad Ka emission line, the latter as required by the si- 
multaneous RXTB-PCA data, and the FeXXV and CaXIX/XX absorption lines. 



Taking pile-up into account (Section 2.2), the nondip Chan- 
dra spectra alone can already be described quite well by a 
weakly absorbed, relatively flat power-law spectrum with a pho- 
ton index Ti = 1.51 ± 0.01 (Table 2). This result is consistent 
with the fact that the break energy of the broadband broken 
power-law spectrum is found at £break = 9keV, i.e., the Chan- 
dra data are virtually entirely in the regime of the (steeper) 
photon index F\. The onset of the exponential cutoff (with 
folding energy Zsfoid) is at E cut = 24keV and thus also well 
above the spectral range of Chandra. As it is known that there 
are cross-calibration uncertainties between Chandra and RXTE 
(Kirsch et al. 2005), we use constant factors c, for the relative 
flux calibration of every spectrum and also separate parame- 
ters Ti(HETGS) and Ti(PCA), for which we find similar values 
within the joint model (see Table 2). In order to describe a weak 
soft excess, we add a thermal disk component, which accounts 
for ~9% of the unabsorbed 0.5-10 keV flux. The disk has a 
color temperature of kT C0 \ = 0.25 keV; similar to that found by 
Makishima et al. (2008) and Balucihska-Church et al. (1995), 
who described the soft excess with a kT = 0.13±0.02keV 
blackbody only, but the temperature T co \ may be too high by 
a factor / > 1.7 (Shimura & Takahara 1995). The norm param- 
eter of the diskbb model is 



'disk 



n . i in ( ^col/km V 

= diskbb . norm = — ; 

V^/lOkpc J 



(6) 



where d is the distance and 8 is the inclination of the disk, which 
can deviate from the orbital inclination / w 35° (Herrero et al. 
1995) as the disk may be precessing with a tilt S 37° 
(Brocksopp et al. 1999). R co i is related to the inner radius of the 
disk, R in = i]g(i)f 2 -R col (with 77 w 0.6-0.7 and g(i) w 0.7-0.8; 
Merloni et al. 2000). In spite of these uncertainties and the 
large statistical error of Adi s k (more than 50 %), the inner disk 
radius can be estimated to 1.5 Rs < R m < lORs if a distance 
d f=a 2.5 kpc (Ninkov etal. 1987a) and a Schwarzschild radius 
Rs i=3 30km (Herrero et al. 1995) are assumed. Thus, the disk is 
consistent with extending close to the innermost stable circular 
orbit (ISCO). 

The good S/N afforded by the RXTE-PCA data clearly re- 
veals an iron fluorescence Ka line. While the instrumental re- 
3.3fcponse of the proportional counters does not allow for a res- 
olution of the line profile details, i.e., whether it is narrow or 
relativistically broadened, the Qiandra-HEG spectra (Figure 6) 
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do resolve a strong narrow component at 6.4 keV. Nevertheless, 
since the integrated flux of the Chandra measured line is insuf- 
ficient to account for all of the PCA residuals in this region, 
we include an additional broad feature to our modeling, which 
is also compatible with the Chandra spectrum. Given the rel- 
atively low S/N at these energies, we do not model the broad 
iron line with a proper physical model such as a relativisti- 
cally broadened line, but use a Gaussian with its energy fixed at 
6.4 ke V, as the latter is hardly constrained by the data. These re- 
sults illustrate the synergy of the simultaneous observation with 
complementary instruments, as the combination of narrow and 
broad line could only be revealed by the analysis with a joint 
model. 

Our global model for the continuum spectrum now enables us 
to address the features of the high-resolution Chandra spectra, 
which is the topic in the remainder of this section. 

3.3. Neutral Absorption 

Absorbing columns can be measured most accurately from 
the discrete edges in high-resolution spectra at the ionization 
thresholds. We detect the most prominent L-shell absorption 
edge of iron and the K-shell absorption edges of oxygen and 
neon. Juett et al. (2004, 2006a) have inferred the fine struc- 
ture at those edges from earlier Chandra -HETGS observations 
of CygX-1 and other bright X-ray binaries: the Fe L2 and L3 
edges, due to the ionization of a 2p\/ 2 or a 2p 3/ / 2 electron, re- 
spectively, are separately detectable at 17.2 A and 17.5 A. The 
O K-edge at 22.8 A is accompanied by (ls^2p) Ka and higher 
(ls^ np) resonance absorption lines. The Ka line occurs at 
23.5 A for neutral O I and at lower wavelengths for ionized oxy- 
gen. In the case of neon, neutral atoms have closed L-shells, 
such that Nel only shows a (ls^3p) Kf3 absorption line close 
to the K-edge at 14.3 A, while ionized neon also shows Ka ab- 
sorption lines, e.g., Nell at 14.6 A and Nem at 14.5 A. Im- 
proved modeling of the neutral absorption that takes these fea- 
tures into account has recently been included in the photoab- 
sorption model tbnew 12 (Juett et al. 2006b), an extension of 
the commonly used tbvarabs model (Wilms et al. 2000). 

As part of the spectral model for the whole continuum, the 
tbnew model can be used to describe the absorption edges 
detected with the Chandra observation of CygX-1 discussed 
in this paper. Figure 7 shows the Fe L-edges requiring a 
blueshift by AA/A • c = (540 ± 230) km s" 1 of the tbnew 
model, which relies on the cross sections of metallic iron mea- 
sured by Kortright & Kim (2000). Unlike Schulz et al. (2002), 
Miller et al. (2005) and Juett et al. (2006a) have also found that 
the Fe L3 edge requires a small shift; their mean position 
of maximum optical depth is AFeL, = 17.498 A, but our value 
A F eL 3 = (17.469 ±0.014) A is still lower. The shift could be 
caused by the Doppler effect due to a moving absorber, by a 
modified ionization threshold due to chemical bonds, or ioniza- 
tion of the iron atoms (van Aken & Liebscher 2002). In an anal- 
ysis of the Cyg X-l high/soft and low/hard state, focused on this 
spectral region, J. Lee et al. (2008, in preparation) find that the 
Fe L-edges here can likely be modeled by a heterogeneous com- 
bination of gas and condensed matter of iron in combination 
with oxygen local to the source environment. If, as suggested 
by these authors, the magnitude of the shift is due to molecules 
and/or dust, this shift is one identifying signature of the compo- 
sition and charge state of the condensed state material. Such 
direct Chandra X-ray spectroscopic detection of dust via its 

12 The code for tbnew is available online at http://pulsar.sternwarte.uni- 
erlangen.de/wilms/research/tbabs/ for beta testing and will soon be released. 
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Wavelength [AJ 

Figure 7. Spectral region around the Fe L-edges. The top panel shows the 
model flux in units of 1CT 5 ph. s _I cm -2 A -1 (solid line). The dashed lines de- 
scribe the unshifted model, see the text. The second panel shows the count rate 
of all the four HETGS spectra (MEG± 1 , HEG± 1 ), which have been combined 
with a resolution of 10 mA, and the folded model. The bottom panel displays 
the residuals \ = (data-model)/error. 
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Figure 8. Same as in Figure 7, but for the spectral region around the O K-edge. 
Both MEGdz 1 spectra have been combined with a resolution of 0. 1 A. 
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Figure 9. Same as in Figure 7, but for the spectral region around the Ne K- 
edge. All spectra have been combined with a resolution of 5 mA. The dotted 
lines show the absorbed continuum model without additional absorption lines. 
The light dashed line shows the model without the absorption of Ne. 
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Table 3 

Column Density N and Abundance A=N / /% of Neutral Absorbers Detected Along the Line of Sight Toward Cyg X-l 



Analysis by 



This work 



ObsID3814 



Schulz et al. (2002) 



ObsID 107 



Juett et al. (2004, 2006a) 



ObsID 107 



ObsID 3407 



ObsID 3724 







(Zfl 


5-10kev = 0.7X10 37 


ergs ') 


(^0.5-lOkev = 


l.exlO^ergs" 1 ) 


(3.1xl0 37 ergs-') 


(soft state) 


Element 


12 + logA 1 ™ 

o '^element 


1 'clcniclll/ "rt 


element / element 


(10 17 cnT 2 ) 


^element'' 
(10 17 cm" 2 ) 


Welemem'' 
(10 17 cnT 2 ) 


J 'clement 

(10 17 cm~ 2 ) 


Welemem'' 

(10 17 cm~ 2 ) 


O 


8.69 


15.4±0.3 


1.45 ±0.01 


38.1^ 


39.2 ±2.3 


63 ±14 


24±3 


29±3 


Ne 
Na 


7.94 
6.16 


3.91 ±0.18 
1.26±0.15 


2.07+S^ 

-0.03 

40±2 


9.6 ±0.2 
3.1±0.2 


9.43 ±0.32 


7.1 ± 1.0 


7.4^ 7 

' -0.3 


8.6 ±0.7 


Mg 


7.40 


2.16±0.13 


4.0±0.1 


5.3±0.2 


3.7 ± 1.3 








Al 


6.33 


<0.13 


< 1.8 


<0.2 










Si 


7.27 


1.45 ±0.17 


3.8±0.3 


3.8±0.3 


2.3 ± 1.8 








S 


7.09 


0.6 ±0.2 


2 4+0.5 


1 6+ 03 










Ar 


6.41 


<0.06 


< 1.0 


<0.1 










Ca 


6.20 


<0.05 


< 1.3 


<0.1 










Cr 
Fe 


5.51 
7.43 


<0.18 
1 


1 3+ 27 
1.75 ±0.03 


02+ 05 

u ' uz -0.02 

2.52 ±0.04 


1.5 ±0.1 


1 47+O.25 


1 06+O 07 

1-UD -0.11 


0.96 ±0.09 


Notes. 



" We use elemental abundances in the interstellar medium, A'^Jj according to Wilms et al. (2000), i.e., xspec_abund ( "wilm" ) ; in ISIS. 
* The relative elemental abundance, A c i cracnt /A'^ 1 i: , is a parameter of the tbnew model. 
C ^element i s the product of A e i emellt and Nh (see Table 2). 

'' Schulz et al. (2002) use the cross sections from Verner et al. (1993), and Kortright & Kim (2000) for Fe. 

" Juett et al. (2004, 2006a) use the cross section of Gorczyca & McLaughlin (2000) for O, Gorczyca (2000) for Ne, and Kortright & Kim (2000) for Fe. 



associated edge structure was first suggested for observations 
of the Fe L-edge in the active galactic nucleus MCG-6-30-15 
(Lee et al. 2001) and for the observations of the Si K-edge in 
the microquasar GRS 1915+105 (Leeetal. 2002). The latter 
study associated the observed Si K-edge structure with Si02, 
although the origin - source environment or Chandra CCD gate 
structure - was unclear in this case. 

The blueshift of (100± 130)kms _I which has been deter- 
mined for all other absorption edges is consistent with zero. The 
O K-edge can only be seen in the MEG spectra after heavy re- 
binning (Figure 8). Nevertheless, the Ka resonance absorption 
line of O I is clearly detected. The region around the Ne edge 
(Figure 9) is dominated by Fe XVIII absorption lines, possibly 
blending with the Ka absorption lines of Ne II and Ne III. No 
other strong edges are clearly visible in the spectrum. The Na 
K-edge (at 1 1.5 A; Verner & Yakovlev 1995) blends with an ab- 
sorption line due to the 2s^3p excitation of Fe XXII. The Mg 
K-edge (at 9.5 A) blends with the Ne X Ly S absorption line. 
The Si K-edge (at 6.7 A) is strongly affected by pile-up and 
blends with the Mg XII Ly (3 absorption line. The S K-edge (at 
5.0 A) is relatively weak. Neutral absorption from these ele- 
ments is nevertheless required within the tbnew model. 

The results for the individual abundances, A, = N,/Nh, and re- 
sulting column densities, A 7 ,, are presented in Table 3. The alpha 
process elements O, Ne, Mg, Si, and S are overabundant with 
respect to the interstellar medium (ISM) abundances as summa- 
rized by Wilms et al. (2000) and therefore suggest an origin in 
the system itself. The total column densities are also compared 
with the values obtained by Schulz et al. (2002) and Juett et al. 
(2004, 2006a) from other Chandra observations of Cyg X-l. 



Table 3 includes the corresponding source luminosities if they 
are reported in the literature (see also Section 4.4). The X-ray 
flux was highest during the soft state observation with the ob- 
servation identification (ObsID) 3724. The column densities 
confirm the conjecture of Juett et al. (2004) that a higher (soft) 
X-ray flux ionizes material local to the Cyg X-l system and re- 
duces the neutral abundances. 

The inferred hydrogen column density A^h (see Table 2) is in 
very good agreement with that from ASCA observations dur- 
ing the soft state in 1996, namely N H = (5.3 ± 0.2) x 10 21 cm" 2 
(Dotani et al. 1997), and also with A^h = 6.2x 10 21 cm" 2 obtained 
from two other different Chandra observations (Schulz et al. 
2002; Miller et al. 2002, see also Section 4.4). We note that 
the large column density toward Cyg X-l found by many on- 
line tools, (7.2-7.8) x 10 21 cm" 2 , is obtained from a coarse grid 
- with (0.675 deg) 2 pixel size - of A^h measurements at 21 cm 
(Kalberla et al. 2005), which does not resolve the strong varia- 
tions of A^h in the region around Cyg X-l (Russell et al. 2007). 



3.4. Absorption Lines ofH- and He -Like Ions 
The high-resolution spectra reveal a large number of absorp- 
tion lines of highly ionized ions. The 1.5-15 A range is shown 
in Figure 10 as the ratio between the data and the continuum- 
model. As the line profiles are not fully resolved, we model 
each line / with a Gaussian profile G/. In terms of the contin- 
uum flux model, F cont , the global model reads 



F(X) = e" T(A) ■ F cont (A) = 



1 + ^G/(A) 



Fcont(A) . (7) 
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Figure 10. Chandra (nondip) spectrum of CygX-1, shown as ratio of data and continuum model (Table 2). For visual clarity, the data have been rebinned to a 
common resolution of 10 mA, and all MEG±1 and HEG±1 spectra have been combined. The Gaussian line fits and identifications are shown as well. (The labels 
mark the lines' rest wavelengths.) Iron L-shell transitions (lines of Fe <XXV) are shown in blue in the online journal. 
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Table 4 

Overview on the Detected Lines from H- and He-Like Ions: Theoretical Rest Wavelengths in A 





Transition 


O 


Ne 


Na 


Mg 


Al 


Si 


S 


Ar 


Ca 


Fe 


Ni 




Hydrogen-like (1 electron) 


VIII 


X 


XI 


XII 


XIII 


XIV 


XVI 


XVIII 


XX 


XXVI 


XXVIII 




ls( 2 S 1/2 )- 


+ 2p( 2 P 3/2 , 1/2 ) 


18.97 


12.13 


10.03 


8.42 


7.17 


6.18 


4.73 


3.73 


3.02 




1.53 


Ly/3 


ls( 2 S 1/2 )- 


+ 3p( 2 P 3/2il/2 ) 


16.01 


10.24 


8.46 


7.11 


6.05 


5.22 


3.99 


3.15 


(2.55) 


1.50 


1.29 


Ly7 


ls( 2 S 1/2 )- 


+ 4p( 2 P 3/2 , 1/2 ) 


15.18 


9.71 


8.02 


(6.74) 


5.74 


4.95 


3.78 


2.99 


2.42 


1.43 


1.23 


Ly5 


ls( 2 S 1/2 )- 


-» 5 P( 2p 3/2,l/2) 


14.82 


9.48 


7.83 


6.58 


(5.60) 


(4.83) 


3.70 


(2.92) 


(2.36) 


(1.39) 


(1.20) 


Lye 


ls( 2 S 1/2 )- 


->6p( 2 P 3 / 2il/2 ) 


(14.63) 


9.36 


( 7.73) 


6.50 


(5.53) 


(4.77) 


(3.65) 


(2.88) 


(2.33) 


(1.37) 




LyC 


ls( 2 S 1/2 )- 


^7p( 2 P 3/2-1/2 ) 


(14.52) 


9.29 


( 7.68) 


6.45 


(5.49) 


(4.73) 


(3.62) 


(2.86) 


(2.31) 


(1.36) 




Lyr? 


ls( 2 S 1/2 )- 


-8p( 2 P 3/2 , 1/2 ) 


(14.45) 


9.25 


( 7.64) 


(6.42) 


(5.47) 


(4.71) 


(3.60) 


(2.85) 


(2.30) 


(1.36) 




Lye 


ls( 2 S 1/2 )- 


-» 9p ( 2 P 3 / 2 ,i/ 2 ) 


(14.41) 


( 9.22) 


(7.61) 


(6.40) 


(5.45) 


(4.70) 


(3.59) 


(2.84) 


(2.29) 


(1.35) 




limit 


ls( 2 S 1/2 )- 


-> oo 


(14.23) 


(9.10) 


(7.52) 


6.32 


5.38 


(4.64) 


(3.55) 


(2.80) 


(2.27) 


1.34 


(1.15) 





Transition 


O 


Ne 


Na 


Mg 


Al 


Si 


S 


Ar 


Ca 


Fe 


Ni 




Helium-like (2 electrons) 


VII 


IX 


X 


XI 


XII 


XIII 


XV 


XVII 


XIX 


XXV 


XXVII 


f [em.] 


ls 2 ('S )+- 


ls2s( J Si) 


22.10 


(13.70) 


11.19 


9.31 


7.87 


6.74 


5.10 


(3.99) 


(3.21) 


(1.87) 




i [em.] 


ls 2 ('S )<- 


ls2p( 3 P 1 , 2 ) 


21.80 


(13.55) 


11.08 


9.23 


7.81 


(6.69) 


5.07 


3.97 


3.19 


(1.86) 


1.60 


r = He a 


ls 2 ('S )^ 


ls2p('P[) 


21.60 


13.45 


11.00 


9.17 


7.76 


6.65 


5.04 


3.95 


3.18 


1.85 


1.59 


He/3 


Is 2 ('So)-* 


ls3p('P,) 


18.63 


11.54 


9.43 


7.85 


6.64 


5.68 


4.30 


3.37 


2.71 


1.57 


(1.35) 


He 7 


ls 2 ( 1 S )^ 


ls4p( 1 P 1 ) 


(17.77) 


11.00 


8.98 


7.47 


6.31 


5.40 


4.09 


3.20 


(2.57) 


1.50 


(1.28) 


He 5 


ls 2 ('S )^ 


ls5p( 1 P 1 ) 


(17.40) 


10.77 


8.79 


7.31 


(6.18) 


(5.29) 


4.00 


(3.13) 


(2.51) 


1.46 


(1.25) 


Hee 


ls 2 ('S )^ 


ls6p( 1 Pi) 


(17.20) 


10.64 


( 8.69) 


(7.22) 


(6.10) 


5.22 


3.95 


(3.10) 








HeC 


ls 2 ('S )^ 


ls7p( 1 P 1 ) 


(17.09) 


10.56 


( 8.63) 


(7.17) 


(6.06) 


(5.19) 


(3.92) 










He r\ 


ls 2 ('S ) -> 


IsSpC'Pi) 


(17.01) 


(10.51) 


( 8.59) 


(7.14) 


(6.03) 


(5.16) 


(3.90) 










limit 


ls 2 ('S )^ 


Is oo 


16.77 


10.37 


8.46 


(7.04) 


5.94 


5.09 


(3.85) 


(3.01) 


2.42 


1.40 


(1.20) 



Notes. Lines with (wavelengths in parentheses) are not detected in our Cnandra-HETGS observation of CygX-1, while lines indicated with bold wavelengths are 
clearly detected and those with underlined wavelengths are detected as two components. The wavelengths of the lines are taken from the CXC atomic database 
ATOMDB and the table of Verner et al. (1996), those of the series limits (= K-ionization thresholds) are from Verner & Yakovlev (1995). 



Table 5 

Results from the Detected Absorption Lines from H- and He-Like Ions: Velocity Shifts v = (A — Aq)/Aq ■ c in kms 



H-like 


O 
VIII 


Ne 
X 


Na 
XI 


Mg 
XII 


Al 

XIII 


Si 
XIV 


s 

XVI 


Ar 
XVIII 


Ca 

XX 


Fe 
XXVI 


Ly a 
Ly/3 
Ly 7 
Ly 5 
Lye 


7 i o+1716 

-90+ 187 
yu -185 

-109+ 1109 


-128+H 

77+50 
/z -47 

7+64 
'-64 

_1Q7+46 
2S8+ 2 " 


284+JI 


-34+ 26 

, 7 +27 
i ;_ 73 

3+922 
-1401 


-171+ 88 
l ' 1 -104 

^44+431 


-60+J 2 

7 c+263 

-238 

308+333 


_43+ii4 

_n 47 +393 

-399 

420+536 


379+712 
Jzz -617 

-750+™ 


164 +838 




Ly series 


-41+67 


_Q7+28 


207+104 
zu '-101 


-38+| 8 


-183+]* 


-54+ 33 


-90+ 74 

w -73 


"327^28 









O 


Ne 


Na 


Mg 


Al 


Si 


s 


Ar 


Ca 


Fe 


He-like 


VII 


IX 


X 


XI 


XII 


XIII 


XV 


XVII 


XIX 


XXV 


He a 




-203+ 8 ] 




-71-11 


286+« 7 


-30+ 23 
^-46 


04+127 
°^-138 


_, 42+250 
1 ^ z -233 


1 1 24+ 98 
llz ^-1069 


11D -531 


He/3 


,07+492 
JU '-1020 


-36+f 




-65^ 




1 S+165 
1S -195 




77+440 
' -306 


_74+616 




He 7 




-507+ 106 
Ju -123 




77 +8l 
' '-71 




1 q+320 
1 -234 


702+™ 


65 5+ 283 






He<5 






















Hee 




-93+g 


















He series 


-291 ±(> 1000) 


-158+35 




_ 72 +25 

/z -28 


135^| 


-86+26 




7+8 
z -12 




-21+° 

-379 



Notes. A negative velocity indicates a blue shift, due to the absorbing material moving toward the observer. Rows labeled with "Ly/He series" show the results from 
modeling the complete absorption line series of the corresponding ion at once with a single physical model (Section 3.5). 
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Table 6 

Results from the Detected Absorption Lines from H- and He-Like Ions: Column Densities in 10 16 cm 





O 


Ne 


Na 


Mg 


Al 


Si 


S 


Ar 


Ca Fe 


H-like 


VIII 


X 


XI 


XII 


XIII 


XIV 


XVI 


XVIII 


XX XXVI 


Ly a 




3 5+ 02 
J - -0.3 




5-2=8:1 


1.8±0.5 


7-1^ 




2 2 +L2 


2±2 


Ly/3 


15« 


10.8+J; 2 




8.(tiJ 


3+ 3 

--2 


18« 




7+?, 




Ly 7 




513 








19 ± 13 


18^ 






Ly <5 


1 1+13 






o+18 
S -8 






9Q+34 
zy -29 






Ly series 


36 +18 


42+10 


7 n+0.4 


7.2 ±0.6 


1.4 ±0.5 


10.1 ±0.8 


l5 + _\ 2 


"3 







O 


Ne 


Na 


Mg 


Al 


Si 


S 


Ar 


Ca 


Fe 


He-like 


VII 


IX 


X 


XI 


XII 


XIII 


XV 


XVII 


XIX 


XXV 


He a 




1.5 ±0.2 




'■-"-0.15 


°-5ti 




<-> a+0.4 
z - J -0.5 


1.3 ±0.6 


1 o +0 9 




He/3 


*3 






5.6±0.9 


3-8^ 


n 


16* 


4±4 


11+ 5 

u -7 




He 7 




3.0 ± 1.5 




8±2 






23!!i 


2l!!o 






He series 


0^ 


6 3 +L3 




5 +2 
-i 




12±3 


15-1 


8^ 




146^° 



Notes. The column densities for the single lines have been calculated using Equation (11), assuming that the line is on the linear part of the curve of growth, which 
underestimates the column density for saturated lines. Rows labeled with "Ly/He series" show the results from modeling complete absorption line series (Section 3.5). 



From a Gaussian's centroid wavelength A and the rest wave- 
length Aq of the identified line, the radial velocity 



v = (A-Ao)/A - 



(8) 



of the corresponding absorber can be inferred. With the defini- 
tion in Equation (7), the norm of G/ is just the equivalent width: 



dA = G,(X)d\ (9) 



W\ } i is related to the absorber's column density, as a bound- 
bound transition i — > j (with the rest frequency vq and oscillator 
strength fij) in an absorbing plasma with column density A/,- cre- 
ates the following line profile (see Mihalas 1978, Section9-2): 



e"™ = exp 



fr 



H 



u-u 



m e c Ai/d \47rAfD At'D 



(10) 



Assuming pure radiation damping, 
equals the Einstein coefficient Aj,. 



the damping constant T 
The Doppler broadening 



Az/d = i/ ■ £ /c is given by £ = 2^r/m ion +vf urb , i.e., is due to 
the thermal and turbulent velocities of the plasma. For optically 
thin lines (with t{v) < 1), the equivalent width is independent 
of Az/d, such that the absorbing column density can be inferred 
(Spitzer 1978, Eq. 3-48; Mihalas 1978, §10-3): 



N t = 



mc ■ W\ 



1.13xl0 17 cirT 2 



fj 



5> 



If the lines are, however, saturated, W\ depends on Ai/u as well, 
and one has to construct the full "curve of growth" with several 
lines from a common ground state ; in order to constrain N/ (see, 
e.g., Kotani et al. 2000). 

We have therefore performed a systematic analysis of absorp- 
tion line series: H-like ions are detected by their Is — > np Ly- 
man series and He-like ions by their Is 2 — > lsnp resonance 
absorption series, see Table 4. For those lines that are clearly 
detected and not obviously affected by blends, the measured ve- 
locity shifts (Equation 8) are shown in Table 5. Most of the lines 
are detected at rather low projected velocity (|v| <200kms _1 ). 
Note that the systemic velocity of Cyg X-l/HDE 226868 is 



(-7.0±0.5)kms _I (Gies et al. 2003), and that the radial veloc- 
ity of both the supergiant and the black hole vanishes at or- 
bital phase (/> 01 b = 0, while the radial component of the focused 
stream should be maximal (likely to be up to 720 km s -1 , see 
Section 4. 1). The column densities in Table 6 are calculated us- 
ing Equation (11), assuming that the line is on the linear part 
of the curve of growth. As the strongest lines are, however, 
often saturated, Equation (11) predicts too small column densi- 
ties from them. For weak lines, however, the equivalent width 
is most likely to be overestimated such that the quoted values 
may rather be upper limits. The properties of the lines (Einstein 
A-coefficients and quantum multiplicities, which determine the 
oscillator strength, as well as the rest wavelengths) are taken 
from Verner et al. (1996) and ATOMDB 13 version 1.3.1. 

3.5. Line Series of H-/He-Like and Fe L-Shell Ions 

As an alternative to the curve of growth, we chose to develop 
a model which implements the expected line profiles (Equa- 
tion 10) directly for all transitions of a series from a common 
ground state i. The model contains N,, £o, and the systemic 
shift velocity v (Equation 8) as fit parameters, and thus avoids 
the use of equivalent widths at all. This approach allows for 
a systematic treatment of the iron L-shell transitions as well, 
which are often blended with other lines such that the different 
contributions to a line's equivalent width can hardly be sepa- 
rated when only single Gaussians are used. As an example, 
Figure 1 1 shows the Fe XIX complex between 12.8 A and 14 A. 
Lines from different Fe XIX transitions overlap, and so does the 
strong Ne IX r-line. Furthermore, the absorption features are of- 
ten rather weak and no prominent lines can be fitted, whereas 
the line series model can still be applied. Although a disad- 
vantage of this approach is the larger computational effort, it 
usually allows us to constrain the parameters of a line series 
more tightly. 

The model with physical absorption line series fits the data 
hardly worse than the model with single Gaussian lines: x 2 of 
12812 instead of 12180 before (see Table 2) is obtained. The 
results are presented in the last row for the whole Ly/He series 
of Tables 5 and 6 for the H-/He-like ions, and in Table 7 for the 

13 See http://cxc.harvard.edu/atomdb/. 
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Figure 11. Absorption line series of FeXlX between 12.8 A and 14 A. The 
Ne IX r line which blends with the complex at 13.4-13.6 A is also shown, and 
the expected positions of the Ne IX i and f emission lines are indicated. 

Table 7 

Parameters of the Absorption Line Series for Fe L-Shell Ions 



Ion 


V 


W(ion) 






(kms- 1 ) 


(10 16 cnr 2 ) 


(kms" 1 ) 


Fe XXIV 


21+ 31 
zl -30 


4 q+0.7 
+■^-1.0 




Fe XXIII 


120+g 


9 9+0.4 
z - z -0.2 


410% 


FeXXII 


-27+* 


2 3+°' 2 




FeXXI 


-156+I" 


2 5+°' 2 
"-0.4 


260+J* 


FeXX 


-283+33 


3 0+ - 3 


333+g 


FeXIX 


-116+* 


3 1+ 02 


251+|° 


Fe XVIII 


-60+* 


LS -0.2 


203g 2 


Fe XVII 


-277+|f 


1.4 ±0.3 


78+| 



Fe L-shell ions. The column densities inferred from the series 
model are generally in good agreement with the values derived 
from the single Gaussian fits for the higher transitions of H- 
or He-like ions (Table 6), while the a - sometimes even f3 - 
lines are saturated. These measurements can be used as input 
for wind or photoionization models, although a few columns 
are rather badly constrained if the thermal velocity £o is left as 
a free parameter. As the line profiles are not resolved, there 
is a degeneracy between Doppler-broadened lines and narrow, 
but saturated lines. The notable fact that no strong wavelength 
shifts are observed is, however, independent of this degener- 
acy: almost all line series are consistent with a velocity between 
-200 kms" 1 and +200 kms" 1 . Lower ionization stages of the 
same element are usually seen at higher blueshift, like in the 
sequence Fe xxiii-xxn-xxi-xx. 

3.6. Emission Lines from He-Like Ions 

The transitions between the Is 2 ground state and the ls2s or 
ls2p excited states of He-like ions lead to the triplet of for- 
bidden (f), intercombination (i), and resonance (r) line, see Ta- 
ble 4. These lines provide a density and temperature diagnos- 
tics of an emitting plasma via the ratios R(n e ) := F(f) / F(i) 
and G(T e ) := [F(i) + F(f)] / F(r) of the fluxes F in the r-, i-, 
and f-line (see, e.g., Gabriel & Jordan 1969; Porquet & Dubau 
2000). In this observation, the dipole-allowed resonance transi- 
tions are seen as absorption lines, as an absorbing plasma is 




9.15 9.2 9.25 9.3 9.35 



Wavelength [A] 

Figure 12. Mg XI triplet of resonance (r), intercombination (i), and forbidden 
(f) lines in the Ciiandra-HETGS spectrum. The MEG data are shown with 
an offset of —0.5. While the ion/line labels indicate the rest wavelengths, the 
abbreviations rn, in, and fn label the actual line positions (see Table 8). 



Table 8 

Parameters of the Lines of the Mg XI Triplet 



Line 


Component 


A 

(A) 


V 

(kms" 1 ) 


W X 
(mA) 


MgXI r 
Mg XI r? 


rl 

(r2) 


-0.002 

q iqn+0004 
yl -0.003 


-70+35 

(700$°) 


-8.6+™ 
-3.5 ±0.9 


MgXI i 
MgXI i 


il 
i2 


9.232 ±0.003 

9 261+ 002 
' " -0.003 


968!??o 


+ 1 5+ u 
+ 1 "'-0.6 

+1 7+ M 
+ 1 -'-0.9 


MgXI f 
Mg XI f 


fl 

f2 


q ii 1+0.004 
yJ1 -0.001 

q 198+0.009 
" -0.003 


-110+ 130 
45 4+293 


+2 5+°- 9 
+ Z - J -1.2 

+ 1 3+°- 9 



Notes. The shift velocity v and equivalent width W\ are given by Equations (8) 
and (9). Note that the (r2) component is likely to be caused by absorption from 
blueshifted transitions of Fe XXI, Fe XX, and Ne X. 

detected in front of the X-ray source. For the same reason, 
we cannot use the Fe L-shell density diagnostics (Mauche et al. 
2005) as, e.g., the Fe XXII emission lines at 11.77 A and 
1 1.92 A used by Mauche et al. (2003) are both seen in absorp- 
tion. But we can still use the detected He-i and -f emission 
lines to estimate the density via the /?-ratio, noting the caveat 
that the densities are systematically overestimated in the pres- 
ence of an external UV radiation field, as photoexcitation of 
the ls2s ( 3 Si) — > ls2p ( 3 P) transition depopulates the upper 
level of the f-line in favor of the i-line and leads to a lower 
tf-ratio (see, e.g., Mewe & Schrijver 1978; Kahnetal. 2001; 
Wojdowski et al. 2008). 

The i- and f-lines of the Mg XI triplet are seen as two distinct 
components each - one with almost no shift, which is consistent 
with the r-absorption line, and the other one at a redshift of 400- 
1000km s" 1 (see Figure 12 and Table 8). Given these two emis- 
sion components and their Doppler shifts, one could be tempted 
to identify the absorption feature at 9.19 A with a second red- 
shifted Mg XI r line, but our model for the complete absorption 
line series (Section 3.5) predicts that the blueshifted 2s — ► 4p 
transitions from the ground states of Fe XXI and Fe XX, as well 
as the Ne X Is — > lOp transition, account for most of the absorp- 
tion seen at 9.19 A. Table 9 shows the R -ratios obtained for the 
two pairs of lines as well as the corresponding densities accord- 
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Table 9 

Electron Densities Corresponding to R = F(f)/F(i) Ratios 







Model Calculations by 
Porquet & Dubau (2000)" 


Component 


Measurement 


R(Mg XI) 


(10 l2 cnr 3 ) 


First pair 




2.3 


1.0... 4.5 


of i and f lines 


Ri(MgXl)=1.6^ | 


1.6 


3.5... 10 


(unshifted) 




1.0 


8 ... 22 


Second pair 




1.3 


5... 15 


of i and f lines 


R 2 (MgXl) = 0.8±0.5 | 


0.8 


12... 28 


(redshifted) 




0.3 


43... 91 



Note. " A temperature between 0.3 and 8 MK is assumed 
and UV-photoexcitation is not considered. 



Table 10 

Solutions to Equation (15) for the Model of Equation (13) 



nn(x) 


d(x) 


fo 




.V 


T 1 f\{\ I't 

2100 -fix) 


(1U cm ) 




- I'o/Voo 




„ / n 

— V 1 1\± 


tor p - pmax 


440 


200 


0.005 


< oo 


1 


11 


110 


50 


0.011 


< 1 


< 1.01 


<41 


110 


50 


0.011 


< 2 


< 1.08 


< 36 


110 


50 


0.011 


< 3 


< 1.18 


< 30 


22 


10 


0.011 


< 1 


< 1.08 


< 180 


22 


10 


0.011 


< 2 


< 1.29 


< 127 


22 


10 


0.011 


< 3 


< 1.48 


< 95 


4.4 


2 


0.011 


< 1 


< 1.36 


< 570 


4.4 


2 


0.011 


< 2 


< 1.69 


< 368 


4.4 


2 


0.011 


< 3 


< 1.97 


< 271 



Note. The binary separation a = 41 Rq corresponds to x a = a/R* = 2.4. 



ing to the calculations of Porquet & Dubau (2000, Fig. 8e), ne- 
glecting the influence of the UV radiation of HDE 226868 on the 
R-ratio. The unshifted lines would then be caused by a plasma 
with an electron density « e .i ~ 5 x 10 12 cm~ 3 ; the redshifted pair 
of lines seems to stem from another plasma component with 
"e,2 3> «e.i- A more detailed discussion of these densities is 
presented in Section 4. 1 . 

The S/N of the spectrum above 20 A is not high enough to 
describe the O VII triplet. The Ne IX triplet blends with sev- 
eral Fexix lines (Figure 11). Nax i- and f-lines are likely to 
be present with several components, but those cannot be distin- 
guished clearly. Similar as for Mg, there are also two i-lines 
of Alxil at shifts of -477!^] km s" 1 and +333^f km s" 1 , re- 
spectively, but no f-line is detected due to a blend with the 
strong Mg XI He f3 absorption line and especially an Fe XXII 
ls 2 2s2p(2s^5p) absorption line at 7.87 A. Similarly, the Si XIII 
i-line is not detected as it overlaps with absorption features 
probably due to nearly neutral Si K-edge structures. (Further- 
more, the flux of the Si XIII f-line might be underestimated as it 
blends with the Mg XII Ly 7 line.) The i- and f-lines of He-like 
sulfur are detected with R(S XV) = 0.61+qIq, but S XV has not 
been modeled by Porquet & Dubau (2000). Below 4 A, the S/N 
below is again not high enough to resolve further triplets from 
heavier elements such as argon, calcium, or iron. 

4. DISCUSSION AND CONCLUSIONS 

In the following, we discuss our results and derive constraints 
on the stellar wind in the accretion region. 

4. 1 . Velocity and Density of the Stellar Wind 

While a spherically symmetric model for the stellar wind 
in the HDE226868/CygX-l system can be excluded by ob- 
servations (see, e.g., Gies & Bolton 1986b; Gies et al. 2003; 
Miller et al. 2005; Gies et al. 2008), a symmetric velocity law 

v(r) = voo • /(//#*) (12) 

is usually assumed to obtain a first estimate of the particle den- 
sity in the wind. The fraction / of the terminal velocity Vqo is 
often parameterized by (Lamers & Leitherer 1993, Equation 3) 

/(*) = f + (l-/ )-(l-l/x) /J (forx>l) (13) 

with fo := Vrj/Voo, where vn is the velocity at the base of the 
wind (x = 1). The simple model for the radiatively driven 
wind of isolated stars by Castor et al. (1975) is obtained for 



(3 = 1/2. The photoionization of the wind, however, sup- 
presses its acceleration (Blondin 1994), such that a smaller 
/ (e.g., a larger /? within the same model) is required in the 
Stromgren zone. Gies & Bolton (1986b) have explained the 
orbital variation of the 4686 A Hell emission line profile of 
CygX-1 with a similar model for the focused wind, where 
Vqo, R*, an d depend on the angle 9 from the binary axis. 
(3 was interpolated between 1 .60 and 1 .05 for 9 between and 
20°. The value (3(9 = 20°) = 1.05 is, however, often used for 
a spherically symmetrical wind as well (Lachowicz et al. 2006; 
Szostek & Zdziarski 2007; Poutanen et al. 2008). Vrtilek et al. 
(2008) use fo = 0.01 to avoid numerical singularities and fit 
the (relatively low) value (3 ~ 0.75 to their models for UV 
lines, vo is likely to be of the order of the thermal veloc- 
ity of H atoms, which is (2kT jm^l 1 = 23kms _1 and corre- 
sponds to fo = 0.011 for HDE226868's effective temperature 
r eff = 32 000 K and = 2100kms _1 (Herrero et al. 1995). 

With the mean molecular weight per H atom, fj, sw 1.4, 14 
the continuity equation = ^imunuir) ■ Anr 2 v{r) gives the 
following estimate for the hydrogen density profile: 

«hW = M * /( £ mH) -d(r/R+) with d(x) = X ~ 2 /f(x) (14) 

Using the parameters of HDE 226868 (M* = 3 x 10 _6 M©yr\ 
R+ = ITRq; Herrero et al. 1995; also summarized by 
Nowak et al. 1999, Table 1), Equation (14) predicts 

n H (r) = 2.2xl0 10 cm" 3 • d^r/R*) . (15) 

Equation (14) can only be solved within the model of Equa- 
tion (13) -or any other model for / in Equation (12) with 
/«>/o-fbr 

d < fo 1 and 1 < x < y/l/(f d) . (16) 

For fo ~ 0.01, the value d w 190 - which would be required to 
explain the density 15 «h.i ~ 4.2 x 10 12 cm" 3 obtained from the 
unshifted MgXl triplet (Section 3.6) - can never be reached 
within our model for the continuous spherical wind. This result 
shows that the density is likely to be overestimated by an R-mtio 
analysis which ignores the strong UV-flux of the 09.7 star. Ta- 
ble 10 lists therefore some solutions to Equation (15) for much 
lower densities as well. Due to the additional constraint that 
the radial velocity is less than 100 km s" 1 (Table 8), we suggest 

14 The helium abundance per H atom is ft: 10% (Wilms et al. 2000). 

15 A plasma of fully ionized H and He contains ssl.2 electrons per H atom. 
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that the first emission component stems from close to the stel- 
lar surface. Although the simple wind model of Equation (13) 
may not be appropriate in this region, the results of Table 10 are 
rather insensitive to the assumptions on the velocity law, as a 
wide range of /3 values was considered, but mostly depend on 
the wind's initial velocity vo, for which we have used a reason- 
able estimate. Note that fo and d also depend on Vqo, but their 
product, which is important in Equation (16), does not. 

In spite of the systematic errors of the (absolute) density anal- 
ysis with the R-rsXio, we infer that the second plasma com- 
ponent is much denser relative to the first one. As it is seen 
at a larger redshift of 400- 1 000 km s" 1 , we favor its identifi- 
cation with the focused wind. The two emission components 
could also be caused coincidentally by dense clumps in the stel- 
lar wind (which are common for O-stars; e.g., Oskinova et al. 
2007; Lepine & Moffat 2008), but the interpretation as a slow 
base of the wind close to the stellar surface and a focused wind 
between the accreting black hole and its donor star provides a 
consistent description of both emission components: an undis- 
turbed wind (with and fo as above, and (3 = 1 .05 — > 1 .6) 
would reach a velocity of v(a) « 900 <— 1200 km s" 1 in the dis- 
tance of the black hole. The focused wind in the orbital plane 
would then be detected with a projected velocity v(a)sin; = 
500 *— 700 km s" 1 at (j> m \, = 0. While photoionization reduces 
the efficiency of acceleration for the spherical wind, the denser 
focused wind is less strongly affected due to self-shielding and 
reaches the expected velocity. It is an observational fact that 
the focused wind has another ionization structure: optical emis- 
sion lines (Ha, He II A4686) from ions which only exist at a 
low ionization parameter have been observed in the focused 
stream (see, e.g., Gies & Bolton 1986b; Ninkov et al. 1987b; 
Gies et al. 2003). For the spherical wind, however, Gies et al. 
(2008) have conjectured that the part between Cyg X-l and the 
donor star might be completely ionized in the soft state. Our 
observations show that the situation in the hard state may be 
similar. 

4.2. Modeling of the Photoionization Zone 

We use the photoionization code XSTAR 2.11n7b 16 
(Kallman & Bautista 2001) to model the photoionization zone. 
As the latter is quite complex (due to the inhomogeneous wind 
density, which is strongly entangled with the X-ray flux), we do 
not claim to describe the photoionized wind self-consistently, 
but only want to derive a first approximation. For an optically 
thin plasma, the relative population of a given atom's ions is 
merely a function of the ionization parameters 



£to = — 2 

«h r l 



"13 ^2 



erg cm s 



(17) 



where L37 is the ionizing source luminosity above 13.6 eV 



in 10 37 erg s , 



«i3 is the hydrogen density in 10 cm and 



rn is the distance from the source in 10 12 cm. 

For Cyg X-l as observed in the hard state, extrapolation of 
the unabsorbed model obtained by our analysis (Table 2) gives 
L37 « 3.5. Although there are obviously strong variations in 
the density, the XSTAR calculation has been performed with 
constant n^ = 0.01, which is the average of nu(r) from Equa- 
tion (15) for R+ < r < a. Figure 13 shows the resulting popu- 
lation of Fe ions. The ionization parameters at which the pop- 
ulation distributions of some ions of interest peak, as well as 
the corresponding full width at half maximum (FWHM), are 
presented in Table 1 1 . We have also included Si IV, N V, and 

16 See http://heasarc.gsfc.nasa.gov/lheasoft/xstar/xstar.html. 



K-shell ions 



L-shell ions 




3 2.5 2 

i°gi; 

Figure 13. Relative population of the ionization stages of iron as a function of 
the ionization parameter £,=L/ (nu r 2 ), as calculated with XSTAR for L37 = 3.5 
and «i3 = 0.01. 

Table 11 

Ionization Parameters for Peak Ion Populations 



Ion 


log? 


FWHM(logO 


rn= (L 37 /« 13 /10 lo ^) 1/2 


FeXXVI 


2.91 


0.27 


0.6-0.8 


FeXXV 


2.75 


0.32 


0.7-1.0 


Fe XXIV 


2.65 


0.30 


0.8-1.1 


FeXXIII 


2.58 


0.28 


0.8-1.1 


Fe XXII 


2.51 


0.31 


0.9-1.3 


FeXXI 


2.44 


0.36 


0.9-1.4 


FeXX 


2.32 


0.44 


1.0-1.7 


FeXIX 


2.12 


0.40 


1.2-1.9 


FeXVIII 


2.01 


0.23 


1.6-2.0 


FeXVII 


1.97 


0.14 


1.8-2.1 


CaXX 


2.69 


0.56 


0.7-1.2 


CaXIX 


2.39 


0.54 


0.9-1.6 


ArXVIII 


2.55 


0.65 


0.7-1.5 


ArXVII 


2.22 


0.50 


1.0-1.8 


SXVI 


2.38 


0.69 


0.8-1.7 


SXV 


2.08 


0.43 


1.2-2.0 


Si XIV 


2.21 


0.65 


0.9-2.0 


Si XIII 


1.96 


0.39 


1.5-2.4 


SilV 


1.09 


0.18 


4.7-5.7 


MgXII 


2.06 


0.58 


1.2-2.3 


MgXI 


1.82 


0.31 


1.8-2.6 


NeX 


1.94 


0.40 


1.6-2.5 


Ne IX 


1.74 


0.24 


2.1-2.8 


O VIII 


1.80 


0.26 


1.9-2.6 


O VII 


1.62 


0.22 


2.5-3.2 


Nv 


1.51 


0.01 


3.25-3.30 


Civ 


1.48 


0.19 


3.3^1.1 



Notes. We list the range in ionization parameter and distance for the maximum 
ion population (Figure 13) obtained in an XSTAR simulation with L37 = 3.5 
and no = 0.01. Note that the binary separation in units of 10 cm is a\2 = 2.9. 



ClV for a comparison with the work of Vrtileketal. (2008) 
and Gies et al. (2008). The H- and He-like ions detected with 
Chandra only appear at considerable distance from the X-ray 
source. Taking into account that the photons emanating from 
Cyg X-l first propagate through a lower density wind and that 
the (eventually stalled) wind of higher density is only reached 
in the vicinity of the star, the actual distances will be even larger 
than the ryi values quoted in Table 11. 
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Figure 14. Gray-scale image of v ra a .bh(T) of Equation (18), i.e., projected wind 
velocity against the black hole. A spherically symmetric velocity law given by 
Equations (12) and (13) with Voo = 2100kms~ 1 , f a = 0.01, and /3 = 1.05 is 
assumed. The star and the black hole are shown as filled black circles; the size 
of the latter is the accretion radius 2GMbh / v(a) 2 . On the black circle passing 
through the center of the star and through the black hole iVad.BHW is since 
ot(r) = 90° . Positive velocities (+) are seen only within this circle; a lighter gray 
means a larger redshift. Negative velocities (— ) can likewise only occur outside 
of this circle; a lighter gray here means a larger blueshift. The dashed gray lines 
of constant v ra d.BH are shown from —1400 km s to +1000 km s in steps of 
200kms~'. The two highlighted sectors contain all observable lines of sight 
toward the black hole (after rotation in this plane) for the inclination / = 35°. 
The labels show the corresponding orbital phases. 



4.3. Origin of Redshifted X-Ray Absorption Lines 

Additional constraints on the accretion flow can be derived 
by considering the Doppler shifts observed in absorption lines 
(Sections 3.4 and 3.5). Models for the wind velocity like that of 
Equations (12) and (13) predict a velocity v(r) at the position r 
in the stellar wind. But only the projection of v against the black 
hole, v ra d,BH, can eventually be observed as radial velocity in an 
X-ray absorption line. With the angle a(r) between v(r) and the 
direction from r toward the black hole, v ra dBH('") is 



Vrad,BHW = COS Ol{t) ■ |v(/) 



(18) 



Assuming a velocity field with radial symmetry with respect 
to the star, a(r) is just the angle at r between the star and the 
black hole. For example, a is 90° (and v ra d,BH is therefore 0) on 
the sphere containing the center of the star and the black hole 
diametrically opposed. Redshifted absorption lines can only be 
observed from wind material inside this sphere, while the part 
of the wind outside of it is always seen at a blueshift - a fact 
which is independent of the assumed velocity field as long as it 
is directed radially away from the star. 

For a spherically symmetrical wind model, any line of sight 
can be rotated to an equivalent one in a half-plane limited by 
the binary axis. Only a sector with half-opening angle i can be 
observed unless the system's inclination is i = 90°. Figure 14 
shows the projected velocity v m & m ( r ) f° r a simple wind model 
and two sectors containing all observable lines of sight toward 
Cyg X-l for an inclination of i = 35°. 

We now investigate the region where the projected wind ve- 
locity (Figure 14) is compatible with the observed Doppler 
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Figure 15. Orbital phase and mean ASM count rate of the previously reported 
Chandra observations of Cyg X-l (ObsID 3814 is presented in this paper). 

shifts (Tables 5 and 7). We are confident that this method 
allows for the identification of the absorption regions, as the 
low observed velocities are always found close to the a = 90° 
sphere - independent of the wind model. For most of the in- 
vestigated ions (Nex, Mgxii, Sixiv and xm, Sxv, Arxvn, 
Fexix and XVIII) the inferred distance from the black hole 
agrees with the predictions of Table 11, e.g., the projected ve- 
locity -117kms _1 < I'Nex < -69kms _1 measured for Nex is 
(during 0.93 < </> or b < 1) obtained at r n = 1.78 ±0.07. For 
many other ions, both results are still very similar: e.g., Ne IX 
with -214kms _1 < v'Neix < -123kms _1 is expected at r\2 = 
1.95 ±0.07 from the v rac j.BH model and at r\% = 2.1-2.8 from 
the photoionization model. Only for the highly ionized iron 
lines, the small distances are - as already anticipated in Sec- 
tion 4.2 - underestimated by the XSTAR simulation run with 
constant average density, which overestimates the wind den- 
sity close to the black hole. For example, the velocity range 
-9kms _1 < Vpexxrv < 52kms _1 measured for Fexxiv corre- 
sponds to r\2 = 1.51 ±0.07, while the population of this ion 
peaks at ryi = 0.8-1.1 within the model presented in Table 11. 

4.4. Previous Chandra Observations 

In previous Chandra -HETGS observations of Cyg X-l, the 
stellar wind was seen under different viewing angles, or the X- 
ray source was in other spectral states (see Figure 15), which 
probably means that the properties of the wind were also dif- 
ferent. Schulz et al. (2002) have analyzed the 14 ks observation 
ObsID 107 performed in 1999 October (at olb w 0.74), 17 when 
Cyg X-l was in a transitional state. They derive a neutral col- 
umn density of YVh = 6.2 x 10 21 cm -2 from prominent absorp- 
tion edges (see also Table 3) and detect some emission and ab- 
sorption lines with indications of P Cygni profiles. Miller et al. 
(2005) have investigated the focused wind with the 32 ks ob- 
servation ObsID 2415 of Cyg X-l in an intermediate state, per- 
formed in 2001 January (at </> rb ~ 0.77). They report absorp- 
tion and emission lines of H- and He-like resonance lines of 
Ne, Na, Mg, and Si with a mean redshift of salOOkms" 1 , as 
well as some lines of highly ionized Fe and Ni. The column 
density is 6.2 x 10 21 cm" 2 (Miller et al. 2002). Marshall et al. 
(2001) describe Ly a and He a absorption lines, redshifted by 
(450± 150)kms _1 from the 12.6ks observation ObsID 1511 
of Cyg X-l in the hard state, performed in 2001 January (at 

17 Note that Schulz et al. (2002) and Miller et al. (2005) quote an erroneous 
date and thus the wrong orbital phase <j) ^ = 0.93 ^0.74 for this observation. 
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</>orb ~ 0.83). Chang & Cui (2007) report dramatic variability in 
the 30 ks observation ObsID 3407 performed in 2001 October 
(at or b ~ 0.88), when CygX-1 reached its soft state. While a 
large number of absorption lines (mostly redshifted, but not by a 
consistent velocity) is identified in the first part of their observa- 
tion, most of them weaken significantly or cannot be detected 
at all during the second part. The complete ionization of the 
wind due to a sudden density decrement is given as a possible 
explanation. Feng et al. (2003) detect asymmetric absorption 
lines with the 26 ks observation ObsID 3724 of Cyg X-l in the 
soft state, performed in 2002 July (at </> or b = 0). The line cen- 
ters are almost at their rest wavelengths, but the red wings are 
more extended, especially for the transitions of highest ionized 
ions, which they explain by the inflowing focused wind reach- 
ing both the highest redshift and ionization parameter closest to 
the black hole. 

The interpretation of the absorption lines presented in the pre- 
vious section describes our observation consistent with wind 
and photoionization models. It can, however, not be applied 
to all of the other observations; the model of Figure 14 nei- 
ther predicts a conspicuously high redshift at ol t « 0.83 and 
at a higher soft X-ray flux (ObsID 1511), nor a positive red- 
shift at </> rb « 0.77 at a still higher flux (ObsID 2415, see Fig- 
ure 15). Inhomogeneities (e.g., density enhancements in shield- 
ing clumps) and asymmetries of the wind (e.g., due to noniner- 
tial forces in the binary system) may therefore play an important 
role in these cases. 

4.5. Summary 

In this paper, we have presented a Chandra -HETGS observa- 
tion of Cyg X-l during superior conjunction of the black hole, 
which allows us to detect the X-ray absorption signatures of the 
stellar wind of HDE 226868. The light curve near cf> or b = is 
shaped by absorption dips; these are, however, excluded from 
this analysis by selecting nondip times of high count rate only. 
At a flux of ~0.25 Crab, we have to deal with moderate pile- 
up in the grating spectra, for which can be accounted very 
well with the simple_gpile2 model. The continuum of 
both Chandra's soft and RXTE's broadband X-ray spectrum 
has been described consistently by a single model, consisting 
of an empirical broken power-law spectrum with high-energy 
cutoff (which is typical for the hard state) and a subordinated 
disk component with an inner radius close to the ISCO. The 
joint modeling of both spectra reveals the presence of both a 
narrow and a relativistically broadened fluorescence Fe Ka line. 
Chandra has resolved absorption edges of neutral atoms with an 
overabundance of metals, suggesting an origin not only in the 
ISM, but also in the stellar wind of the evolved supergiant. The 
previously suspected antic orrelation of neutral column densi- 
ties and X-ray flux, which is due to photoionization, is con- 
firmed. Absorption lines of highly ionized ions are produced 
where the wind becomes extremely photoionized. For H- and 
He-like ions, Lyman and He-series are detected up to the Ly 
or He e line, which we use to measure the column density and 
velocity of absorbing ions. For the wealth of Fe L-shell transi- 
tions, column densities can best be obtained by directly using 
a physical model for the complete line series of an ion. The 
nondip spectrum shows almost no Doppler shifts, probably in- 
dicating that we have excluded the focused wind by our selec- 
tion of nondip times. We have also detected two plasma com- 
ponents in emission by two pairs of i- and f-lines of He-like 
Mg XI. The first one is roughly at rest and we identify it with 
the base of the spherical wind close to the stellar surface. The 
second plasma component is denser and observed at a redshift 
that is compatible with the focused wind. A simple XSTAR 



simulation indicates that most of the observed ions only exist 
in a distance of the black hole, where the velocity of a spher- 
ical wind, projected against the black hole, is low - which is 
a consistent explanation of the small Doppler shifts observed 
in absorption lines. We review the previously reported Chan- 
dra observations of Cyg X-l and find that not all of them can 
be described in the same picture, as the wind may be affected 
by asymmetries or inhomogeneities. A detailed spectroscopy of 
the absorption dips, which might shed more light on the focused 
wind, will be presented in a subsequent paper. 
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